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0\ , A three-body scattering theory previously proposed by one of the present authors is developed 

^\ • to be applied to the saturated ferromagnetic state in the two-dimensional Hubbard model. The 

0^ ' single-particle Green's function is calculated by taking account of the multiple scattering between 
two electrons and one hole. Several limiting cases are discussed and the relation to the variational 

f^ ' principle is examined. The importance of the three-body correlation is demonstrated in comparison 

^ I with the results of the ladder approximation. A possible phase boundary for the Nagaoka ground 

pL^ ■ state is presented for the square lattice, which improves the previous variational results. 
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I. INTRODUCTION 



. . Electron correlations in narrow bands have attracted much attention after the discovery of high T^ cuprate super- 

conductors. Since a thorough study of electron correlations is quite complicated in actual materials, we often use 
a simpjjfied model to understand the essence of the effects. One of the most widely used models is the Hubbard 
^ ' model,Eiij whose Hamiltonian is given by 

Ti '■ 

C3 ■ i7 = -t ^ aj^fljo- +H.C. + C/^a^-^a^^a^ Oil, (1.1) 

where t > and (i, j) stands for the nearest neighbor pairs. Operator a^o- represents the annihilation of an electron 

I ■ with spin a at site i. 

J>. , In the low density of occupied electrons, the correlation effects on this model are exactly taken jnto account 

T— I ' by considering the multiple scattering between two electrons. Using this low-density theory, KanamoriB studied the 

CN ', ferromagnetism of nickel metal (here holes are in low density) . He found that the effective interaction between electrons 

CN is considerably reduced from the bare value U due to the correlated motion avoiding each other, and that the Stoner 

^^ \ condition for ferromagnetism is strongly modified. In order to extend his theory to higher density, we need to include 

also the effects of the particle-hole multiple scattering, which describes the magnpn excitation on the ferromagnetic 

ground state, and those of the electron-magnon interaction. Hertz and Edwardsu recognized the importance of the 

vertex correction, and determined the electron-magnon vertex by enforcing the Ward- Takahashi identityE However, 

C^ it is impossible to determine the vertex on all points in the momentum-frequency space by the identity alone, and 

C , also it is unclear how to include the particle-particle multiple-scattexkig effect, which becomes important in the low 

I ■ density of electrons. One of the present author has proposed a theoryQ'Ll which takes account of the muJtiolc scattering 

■ O ] between two electrons and one hole with the use of the Faddeev equation of the three-body problem. Oil^ This theory 

5 ■ naturally takes account of the particle-particle and particle-hole multiple-scattering channels on equal footing, and 

thereby it contains the effect interpreted as the vertex correction for the electron-magnon interaction. 

The three-body scattering theory has a clear meaning of the three-body problem on the saturated ferromagnetic 
ground state, although it is applicable to the paramagnetic phase by deriving a similar integral equation to the 
Faddeev equation. In this paper, in order to make clear the physical meaning and the reliability of the approximation, 
5J] ' we confine ourselves to the application to the two-dimensional Hubbard model in the context of the instability of the 
Cu I saturated ferromagnetic ground state. Itphas a nice property inteppolating from the low-densiiy limit to the half-filling 
limit, and satisfies the exact atomic limiilil as previously provedE'cl One of the present authorsQ has already applied the 
three-body scattering theory to the one-dimensional Hubbard model in comparison with the exact solution, having 
calculated the spectral function for the single-particle Green's function. Later Ruckenstein and Schmitt-RinkB have 
examined the same approach in comparison with other theories. However, the actual calculation in two dimensions 
has not been carried out yet. In this paper, we calculate the single-particle Green's function in two dimensions and 
determine the quasi-particle energy. We demonstrate the importance of the three-body correlation by comparison 
with the ladder approximation. We find that a low-energy scale for the quasi-particle energy, which is much smaller 
than its value of the ladder approximation, and that the mass enhancement factor is as large as '^ 10 in the strong 
coupling regime. We discuss several limiting cases for clarifying the approximation made. We also show that the 



present theory satisfies the variational principles This observation helps us to compare the present theory to other 
variational theories. 

As regards the stability of saturated ferromagnetism, Nagaoka considered in his pioneering study the case of one 
hole in an otherwise half- filled band with U = ooOEj He proved that the saturated ferromagnetism is realized as 
ground state on approppatc lattices due to gaining the kinetic enegy of the hole. For two holes on a square lattice 
with U =H?o, Fang et alS3 found that the ground state is a singlet by exact diagonalizations for small clusters. Doucot 
and Wentil found that a long-wave length twisted spin state has lower energy than the saturated ferromagnetic state. 
On the other hand, it has been noticed that finite-size effects are crucial and that closed-shell configurations (numbers 
of holes are 1, 5, 9, • -p^j^vor ferromagnetism whereas open-shell configurations (numbers of holes are 2, 3, At,^ ■ ■ ■) 
tend to destabilize it.E£E3 Quantum Monte Carlo results suggest no ferromagnetic phase for intermediate [/JIj'lJ but 
the cluster sizes used and the simulation temperatures are rather high. Therefore, for a thermodynamic concentration 
of holes, such considerations for a few holes do not serve as proper guides. 

Another route to investigate the stability of the Nagaoka ground state is a variational study for the excitation energy 
with respect to overturning an up-spin electron and placing it on a down-spin band. Of course this approach gives only 
an upper limit for the ferromagnetic boundary. RothtS was the first to discuss in detail the instability of the Nagaoka 
ground state with an approximate scheme equivalent to usipg a variational wave function. It is found superior to the 
Gutzwiller function, since it has more vaciational flexibility.E3 The similar variational wave function was employed by 
Shastry, Krishnamurthy, and Anderson,c2l for discussing the instability of the Nagaoka ground state and the spin- wave 
stiffness. It. was also used to calculate th&,spin-wave spectrum with general momentum in the U = oo limit by Basile 
and ElserEil, and for finite U by Okabe.c3 The softening of the spin-wave energy may set another condition for the 
instability of the Nagaoka state. These calculations leads to an instability for a finite concentration of holes-even in the 
[/ = oo limit. Another interesting variational wave function is constructed by Richmond and Rickayzen,c3 who froze 
the motion of the down-spin electron and solved exactly the states for up-spin electrons under the static potential 
due to the down-spin electron (see Appendix B). This study merely leads to the conclusion that the Nagaoka staiej 
is always stable in the U — oo limit for any concentrations of holes. Recently, Hanish, Uhlig, and Miiller-Hartmancil 
analyzed Roth's wave function and its slight extension on various lattices, and calculated the phase diagram for the 
Nagaoka ground state. The present theory is shown to be equivalent to using a variational wave function similar to 
Roth's but assuming fhc most general form with the states containing one particle-hole pair. It is more flexible than 
that of Hanish et al.f3 and thereby we obtain the phase diagram with a smaller area of the Nagaoka ground state 
than Hanish et al. This may serve as a measure of the accuracy of the present theory. 

In Sec. II, the three-body scattering theory is formulated for the single-particle Green's function, and several limiting 
formulas are derived. The relation to the variational method is also discussed. In Sec. HI, the numerical solution is 
presented for the square lattice. Section IV contains concluding remarks. In the Appendix A, the explicit solution of 
the Faddeev equation is derived in the U — oo limit, and in the Appendix B the method of Richmond and Rickayzen 
is adapted for numerical computations. 

II. THREE-BODY SCATTERING THEORY 

We briefly summarize the three-body scattering theory, and discuss several limiting formulas, which are useful to 
make clear the physical meaning of the approximation made. We also discuss the variational character of the theory. 

A. The Green's function formalism 

We start by assuming the saturated ferromagnetic ground state, where only up-spin states are occupied by electrons. 
No interaction is working between electrons with the same spin. Placing a down-spin electron into the system, we 
study its motion using the single-particle Green's function. 

We make a perturbational expansion, where the Hartree-Fock energy is included into the unperturbed part. Thus 
the Hamiltonian is divided into two parts: 

H^Ho + Hi, (2.1) 

where 

Hf) = ^[eicflkT'^kT + (ek + [/n)aj^j^akx], (2.2) 

k 

^^^J^J2 "L+qi («P2-qT"P2T - («p2-qT"P2T)) Opii- (2-3) 

PiP2q 



Here N is the number of lattice sites and (• ■ •) denotes the average over the ground state. Momentum k is defined 
in the first BriUouin zone, i. e., —n < k^ < tt and — tt < fcy < tt in units of 1/a with a being the nearest neighbor 
distance for the square lattice. The kinetic energy with momentum k is given by ek = —2t{coskx + cos ky) for the 
square lattice. The occupied electron density n is defined by 

nSq.o = J^{Y1 «p-qT"pT)- (2-4) 

p 

The single-particle Green's function is defined by 

G,(k,f) = -*(TK.(t)aL(0)]}, (2.5) 

where T represents the time ordering, and ako-(i) ~ exp(i(i7 — fiNe)t)a]i^e'xjp{~'i{H — fiNe)t) with fi and N^ being 
the chemical potential and the number operator of electrons. On the saturated ferromagnetic ground state, only the 
retarded part does not vanish with the down-spin electron, i. e., 

/•OO 

Gi{Q,z)^~t {T[aQdt)aL{OWdt, z = oj + iij {ri^O+). (2.6) 

Jo 

The unperturbed Green's function is given by 

G}"^(Q,z)-i = z-eQ-f/n + Ai, (2.7) 

with Un being the HF potential. The self-energy is defined by 

G^Q, z)-i = Gf{Ct, z)-i - I]i(Q, z). (2.8) 

We sum up all the Feynman diagrams for the self-energy shown in Fig. |l|. where only three lines of Green's function 
exist in the intermediate states. An important point is that the frequency sum in the intermediate states can be 
explicitly carried out, since the Green's function for the down-spin electron consists only of a retarded part. After the 
frequency sum, we obtain 

Si(Q, z) = (^) { E ^o(Q' ^; pq) + E Go(Q' ^' pq);^Go(Q, z- pqO 

^ ^ pq pqq' 

+ Y. Go(Q, z- pq)Z^Go(Q, z- p'q) + ■••}, (2.9) 

pp'q 

where 

Go(Q, z] pq) = l/(z - eQ+p_q - C/n - Cq + Cp + y). (2.10) 

Hereafter p, p' are restricted inside the Fermi sphere, and q, q' are restricted outside the Fermi sphere. 



We notice that Eq. (2_^) is equivalent to collecting up all the multiple-scattering processes-between two particles 
and one hole. This three-body problem is exactly solvable by using the Faddeev equationnJ We write down the 
Faddeev equation for the scattering matrix T(z); 

T{z)^T^{z)+T2[z), (2.11) 

with 

Here Go(z) is a free propagator l/(z — iJg), ti{z) is the particle-particle scattering matrix, and t2{z) is the particle-hole 
scattering matrix. All these quantities are represented within three-body states. As shown in Fig. 0, we introduce 
three-body states with total momentum Q, which are written as 

|Q;pq)=aJ^+p_q^a]j^apT|i^). (2.13) 



Here \F) is the saturated ferromagnetic ground state. Note that |Q;pq)'s are orthonormal. Then we have 

(Q;pq|^TrlQ;p'q') = Go(Q,z;pq)(5p,p,,5q,qs (2.i4) 

(Q;pq|i.(z)|Q;p'q') ^ ^ ^ _ ^^^(q /p, , + ,^ _ ,) W, (2.15) 



-U 1 

^l + UD2iQ-ci, z - eq + /i) ' 



(Q;pq|i2(^)|Q;p'q') = -^T , , ,,^ ,^ „ ., , , .., ^q.q', (2.16) 



with 



D^ik,z) = ^J2 TJ T^' (2.17) 

A* -^ z — ek-q' — Un — eq> + 2fi 



q 

1 
I?. 



(k,z) = — y ^ . (2.18) 

p/ f 

For solving the Faddeev equation, we introduce the following quantities, 

*Pp(Q,^;p) = 5](Q;pq|Ti(z)|Q;p'q')Go(Q,z;p'q'), (2.19) 

p'q' 

<I>ph(Q,z;q) = ^(Q;pq|T2(z)|Q;p'q')Go(Q,z;p'q'). (2.20) 

p'q' 

It will become clear below that $pp is independent of q and $ph is independent of p. Then the Faddeev equation is 
rewritten as 



^ /o A UDiiQ + p,z + ep-fi) 

*pp(Q,2:;p)- 



l--C/L»i(Q + p,z + ep-/i) 
U 



$ph(Q,z;q) 



l-C/L»i(Q + p,z + ep-/i)iV 
-UD2{Cl-q,z-eq + fi) 



^ Y, Go(Q, z; pq')$ph(Q, z; q'), (2.21) 



q 



l + [/L'2(Q-q,z-eq + yu) 
-U 



1 + UD2iQ - q, z - eq + m) iV -^ 



J-^Go(Q,z;p'q)$pp(Q,z;p'). (2.22) 



An equation containing only $ph(Q, z; q) is obtained by substituting Eq. ( 2.21 ) into Eq. ( 2.22 ). This integral equation 
can be accurately solved for systems of fairly large size, since its variable is only q for fixed Q and z. Using the solution, 
we express the self-energy as 

Si(Q, z) = (^\ Y. Go(Q' ^' pq)[i + *pp(Q' ^; p) + *ph(Q' ^; q)]- (2-23) 

^ ^ pq 

There appear many poles determined from the condition 1 + C/£)2(Q ~ Qj -2 — fq + /^) = in the low-energy region 
for the self-energy. They constitute a continuous spectrum in the single-particle spectra in the N —>■ oo limit. This 
condition is reduced to z = '^Q-q + £q ~ M with Wk being the spin-wave energy in the random phase approximation 
(RPA). The lowest energy of the continuum spectra becomes zero at Q = kj?, since luq + e^^ — /i = 0. For relatively 
small values of \Q\, bound states appear below the continuum. They arc identified as "quasi-particles" , and their 
energy iJq's are determined from 

G(Q, Sq)-i = Sq - [eq - M + Un + E^ (Q, Eq)] = 0. (2.24) 

Now we comment on a widely used approximation called "ladder" approximation, which only considers th e par ticle- 



partic le ladder and particle-hole ladder separately. This is given by neglecting the second terms in Eqs. ( p.21[) and 
(^.22). As a results, the self-energy is given by 



^^^' ^" |7V^l + C/i?2(Q-q,^-eq + A*) iV^l-t/i?i(Q + p,2 + ep-M)j ■ ^^ ^ 

The first term comes from tire particle-hole channel, which is interpreted as the electron-magnon interaction with 
bare vertex U. The second term comes from the particle-particle channel, and its contribution is much smaller than 
the first one in the strong coupling regime (except in the low density). As shown in the next section, the self-energy 
in this approximation is larger in order of magnitude than the present theory. 

B. Limiting cases 

We discuss several limiting cases for clarifying the physical meaning of the present theory. 

1. Low density limit (n -^ Q) 

We expand the following quantities up to the first order of n: 

-UD2iQ-q,z~eq + fi) _ -Un 



l + UD2iQ- q,z - eq + ^) 2; - eq^q - Eq + 2/i' 



(2.26) 



-^ -C/+ ^^- , . (2.27) 



1 + C/£)2(Q - q, z - eq + ^) z-eQ_q-eq + 2^ 



Substituting these relations into Eq. ( ^.22 ), we notice that the first term is order n and so is the second term, if 
$pp(Q,z;p) is order 1. Therefore we can neglect $ph(Q,2:;q) in the lowest order of n, and thereby we obtain the 
lowest-order form of <i>pp(Q, z; p) as 



UDiiQ + p,z + ep-^i) 
l-UDi{Cl + p,z + ep- n)' 



*pp(Q. -; p) - . ;;^^^^':::^'^^^ ^ (2.28) 



Substituting this relation into Eq. (2.23) and using the relation {1/N) ^ Go(Q, z; pq) = Di{Q + p, z + ep — /i), we 
obtain in the lowest order 



^i(Q'-)-^Er 



UDi{Q + p,z + ep-fi) 
Un 



i-c^^Ec 



(2.29) 



N /-^q' 2-eQ_„'-e„'+2M 



Here the chemical potential /i is located at the bottom of the band, so that q' runs oveif,the entire first Brillouin 
zone. This is nothing but the low-density expression Kanamori derived several decades ago.0 Thus the present theory 
covers naturally the exact low-density limit. 

2. Limit to the half-filling (n -^ 1) 



We expand the following quantities up to the first order of 1 — ti: 

C/Z3i(Q + p,z-Kep-^) _ U{l-n) 



1 - UDi{Q, + p, z -I- ep - ^) z- eq+p-q^ + Cp - [/' 



(2.30) 



'^U. .'^"-"' „ . (2.31) 



1 - UDi{Q, + p, z-f ep -/i) z- eq+p-q^ +ep - [/' 

Here qo stands for the mome ntum at the top of the valence band, i. e., qo = (tt, tt) for the square lattice. Substituting 



these relations into Eq. (2.21), we notice that the first term becomes order 1— n and so is the second term, if <i>ph(Q, z; q) 
is order 1. Therefore we can neglect $pp(Q, z; p) in the lowest order of 1 — n, and thereby we obtain the lowest-order 
form of $ph(Q, z; q) as 



*ph(Q,z;q) 



-U—T , 



Eq + p'-qo+Ep'-C/ 



1 + c^^Ep'- 



(2.32) 



P' ^-^Q + p'-qo+S'-'^ 



Substituting this relation into Eq. (2.23) and using a lowest-order relation Go(Q,2:;pq) ~ \/[z 
we obtain up to the first order of 1 — n 



eq+p-q,, + Cp 



U]. 



C/n + E^(Q,z)«L/ 



U{1 



1 + f^^Ep' 



(2.33) 



^P z-Eq + p'-qQ- 



Here p' runs over the entire first Brillouin zone. 

Since available unoccupied states are limited in the up-spin band, the particle-particle multiple-scattering process 
gives rise to merely a higher-order contribution, in contrast to the low -den sity limit. Also, the processes creating more 
particle- hole pairs give rise to higher-oreder contributions. Thus Eq. ( 2.33 ) is exact up to the first order of 1 — n under 
the assumption of the saturated ferromagnetic ground state. But this expression is not so useful, since the saturated 
ferromagnetism is always unstable in this limit by the reason discussed below. 

The second term of Eq. ( 2.33 ) vanishes as n ^ 1, and one may think that the saturated ferromagnetic ground 
state is always stable around n ^ 1. This is not true, since the second term of Eq. (2.33) has a pole at z = wq-q^^, 
which is negative: Sj^(Q, z) ex l/{z — t^Jq-q ) for z ~ i^q-q ■ Here Wk represents the RPA spin-wave energy at n = 1. 
Therefore G|(Q, z)^^ can be zero for z ~ a;q_q^, even if 1 — n is very small. This means that the system is unstable 
through the coupling to the unstable spin-wave excitation. Figure ^ shows cjq for several values of U on the square 
lattice, demonstrating that it is always negative. 



3. Atomic limit (t —> 0, U < oo) 



In the limit of i — > 0, we have 

Go(Q,2;;pq) 



1 



-, Di{q,z) 



1 — n 



D2{(l,z) 



z — Un z — Un z — Un 

Substituting these relations into Eqs. ( 2.21 ) and ( 2.22| ), we solve the Faddeev equation. The solution is given by 

U{\-n) ^ _ , -Un 



(2.34) 



*pp(Q,2:;p) 



$ph(Q,2;q) 



z-U{l~ny ^"'^' '^' z-U{l-ny 
Substituting these relations into Eq. (2.23|), we obtain the self-energy and the Green's function as 



Si(Q,z) = t/ 
Gi(Q,z 



2 n{l-n) 
z-U{l-ny 
1 — n n 



(2.35) 

(2.36) 

(2.37) 



z z-U 

This is the exact form in the atomic limit, first emphasized by HubbardE The fact that tbp three-body scattering 
approximation satisfies the exact atomic limit was first proved by one of the present author.' 



4-. Strong coupling limit (U -^ oo^ 
We expand <I>pp(Q, z\ p) and <I>ph(Q, z; q) in powers of 1/U , 

*PP(Q, z- p) = $(Op)(Q, z; p) + $W(Q, z; p) + 
$ph(Q, z; q) = $("j^)(Q, z; q) + <f>W(Q, z; q) + 



(2.38) 
(2.39) 



Expanding also -Di(k, z) and I?2(k, z) in powers of 1/C/ and substituting them into Eqs. ( 2.21 ) and ( 2.22 ), we derive the 
equations for each order of 1/C/. The derivation of the solutions is straightforward but lengthy. We have summarized 
the derivation in the Appendix A, and here we simply show the result: 



<i>(°)(Q,z;p) = -l, l^ci>a)(Q,,;p) = _^ 



$ 



^^' '^^ J(Q,2;q)' iV^ ph^^'^'^i^ c/ ' 



where 



1 1 



J(Q, z;q) = z-eq + n — ^{eQ+p'-q - v)' 



p' 



F{Q,z)= , ^ , 

AT Z^q J(Q,z;q) 



Substituting these solutions into Eq. ( 2.23 ), and with the help of the expansion, Go(Q,z;pq) 
EQ+p-q — Eq + Cp + /^)/(^'^)^]i wB obtain the self-energy 

I](Q,z) = -C/n-i^(Q,z). 



(2.40) 
(2.41) 

(2.42) 

(2.43) 

-[\l{Un) + {z- 

(2.44) 



The first term cancels out the HF potential. This expression is different from the previously proposed one (Eq. (32) 
in Ref. 8). We verify numerically in the next section that the present expression gives the same quasi-particle energy 
as Roth's variational wave function. 



C. Relation to the variational principle 



We consider a system of A^e up-spin electrons and one down-spin electron with total momentum Q. Let the wave 
function be 



|^(Q))=a|Q)+;^/3pq|Q;pq) 



(2.45) 



Here |Q) = OqAF), and |Q;pq) was defined by Eq. ( 2.13 ), where p and q are restricted inside and outside the 
Fermi sphere in the up-spin band, respectively. Since the expectation value of energy E for this state is given by 
E = {ip\H\'>jj)/{ip\ip), the extremum condition for E is reduced to the following eigenvalue equation: 



Hq-qq + 22 ^Q;pq/3pq = Ea, 
pq 

-ffpq;Qa + 2_^ -ffpq:p'q'/?p'q' = E (3pq, 



(2.46) 



where Hq.q = {Q\H\Q), HQ.p^ = (Q|iy|Q;pq), i:fpq;Q = (Q;pq|iJ|Q), and i?pq;p'q' = (Q; pq|i?|Q; p'q'). Instead 
of solving directly this eigenvalue problem, we calculate the resolvent R{z) — \/{z — H) by using the relation R{z) ■ 
[z — H) = 1. It is given by 



(Q|i?(z)|Q) = 



-ffQ;Q- 2J -f^Q:pq'S'pq;p'q'(^)^p'q';Q 
pq;p'q' 



(2.47) 



where S'pq;p'q'(2:) = (Q;pq|(z — if3)~^|Q; p'q') with H^ being the Hamiltonian operator represented by the compo- 
nents {i?pq;p'q'} within the three-body states. The Faddeev equation is one way of calculating 5'pq;p'q'(z). Therefore, 
(Q|i?(z)|Q) is essentially equivalent to G'(Q,z) except for the origin of energy. Note that the poles of R{z) give the 
energy eigenvalues, and let EiS be such energy eigenvalues. Then the excitation energies of overturning an up-spin 
electron at the Fermi level and placing it on a down-spin band are given by Ei — ^ Cp — /i, which coincides with the 
energies of the poles of G(Q, z). 

Now we compare the present theory with other variational theories. The variational wave function used by Rotht^ 
is given by 



l^^^iQ)) =Y.e^'^^'al il -Y^fit - j)al^a,A a^,^\F), (2.48) 

where kj? is a Fermi momentum, and f{i — j) is tlie variational parameter. The total momentum of this state is 



Q — k^. The wave function used by Shastry et alEB is given by putting /(« — j) — Sij in Eq. (2.48). Roth's wave 
function neglects the correlation between the hole at the Fermi level (described by akj^f) and other excitations. If 
this hole at the Fermi level is disregarded, Roth's wave function consists of three-body states of one particle-hole pair 
in the up-spin band and one electron in the down-spin band. Note that these states are li mited in comparison with 
general three-body states, since two operators are confined on the same site i in Eq. (|2.48|) . This restriction and the 
condition /(O) — 1 are necessary for avoiding the double occupancy in the U ~ oo limit. Therefore, we expect that the 
three-body scattering theory gives the samej-esults as Roth's in the U = oo limit. Actually this is verified numerically 
in the next section. Recently, Hanish el al.c3 analyzed Roth's wave function on various lattices by using the resolvent 
method. They also extended Roth's wave function for finite U by including more doubly occupied states such as 
'Hi'At^jI- Their method has an advantage of being able to treat systems of infinite size, but their wave function is 
still limited in comparison with the general three-body states, since two operators are still confined on the same site 
i. On this point, the three-body scattering theory remains superior to their approach. 

III. CALCULATED RESULTS 

Confining ourselves to finite systems of the square lattice, we calculate the single-particle Green's function. Since 
several ek's are degenerate-for finite systems, shell structures appear on physical quantities, depending on the occupa- 
tion of up-spin electrons JIj'lj The foll owing numerical calculations are carried out on closed-shell configurations. We 



first eliminate $pp(Q, z; p) from Eqs. (2.21) and (2.22) to set up the integral equation for $ph(Q, z; q) only. Then we 
solve the equation on systems of finite size; the matrix of size A^(l — n) x A^(l — n) has to be inverted. 

Figure shows the quasi-particle energy and the lowest boundary of the continuum spectra as a function of Q, for 
6{= 1 — n) = 0.2 and 4t/U = 0.05. The calculation was carried out on a system of A^ = 50 x 50 with a periodic 
boundary condition. The dispersion of the quasi-particle energy Eq is extremely flat in comparison with a free 
dispersion eq. The effective mass m* is defined by Eq = Eq + (l/2m*)Q'^ around Q ~ 0. The mass enhancement 
m* /ra (to is the free mass given by eq) reaches to as large as 12 for this situation. The dispersion curve enters into 
the energy continuum for Qx = Qy > OAn. The quai-particle peak may survive as a resonant peak after entering the 
continuum spectra. With increasing values of |Q|, the boundary of the continuum spectra goes down and eventually 
touches to the x-axis at Q = kp. A small gap seen at Q = k^ in the figure is due to a finite-size effect. 

Figure 0(a) shows the quasi-particle energy at Q = 0, as a function of 4t/C/. The hole density and the system size 
are the same as in Fig. §, (5 = 0.2 and A^ = 50 x 50. The Eq decreases with decreasing values of U. It turns negative 
for 4t/U > 0.085, indicating that the Nagaoka ground state becomes unstable by overturning an up-spin electron at 
the Fermi level and placing it on the bottom of the down-spin band. Figure H(b) shows Eq evaluated in the ladder 



approximation. It is found that the first term in Eq. ( 2.25 ) has a dominant contribution. The absolute values are larger 
in order of magnitude than the values of the three-body scattering theory, suggesting that the ladder approximation 
overestimates the self-energy due to neglecting the vertex correction for the electron-magnon process. Note that the 
three-body scattering theory takes account of the vertex correction, since the particle-hole multiple-scattering process 
is renormalized by the particle-particle multiple-scattering process. 

As already mentioned, the three-body scattering theory is expected to give the same quasi-particle energy as Roth's 
variational wave function in the U = oo limit. We demonstrate this-equivalence in Fig. p), which shows Eq as a 



function of S in comparison with the result of Roth's wave function.tja We have used Eq. (p. 44) for the calculation 



on a system of size A^ = 100 x 100. Both results coincide with each other within the accuracy of numerical errors. 

Now we discuss the phase diagram for the Nagaoka ground state. We determine the phase boundary from the 
condition Eq < 0, which means the instability with respect to overturning an up-spin electron at the Fermi level and 
placing it on the down-spin band. Figure shows the phase diagram thus evaluated in comparison with the results 
of other theories. Note that finite-size effects are not small even for relatively large size A" = 50 x 50. The solid line 
(SW) represents the boundary above which the energy of spin waves becomes negative within the RPA. The solid 
line (RR) represents the boundary determined by the method of Richmond and RickayzenH (see Appendix B). The 
broken line represents the boundary given by Hanish et al. (the curve of "RES3" in Ref . 24) .E£l All the curves have 
the similar slopes around (5 ~ 0, but depart from each other with increasing values of S. The three-body scattering 
theory gives the most severe condition for the phase boundary, the smallest area of the Nagaoka state. The critical 
hole concentration is given by 6c — 0.41, the same as that of Hanish et al. This is reasonable since both theories 
become equivalent to using Roth's variational wave function in the U = cxD-limit. 



IV. CONCLUDING REMARKS 



We have developed a many-body theory which takes account of the three-body multiple scattering, and have 
applied it to the two-dimensional Hubbard model by assuming the saturated ferromagnetic ground state. We have 
solved numerically the Faddccv equation of the three-body problem on systems of finite size, and have calculated 
the single-particle Green's function. We have found that the energy scale of the quasi-particle is smaller in order of 
magnitude than that of the ladder approximation. This difference comes from the lack of the vertex correction to the 
electron-magnon scattering process in the ladder approximation. Recently, a many-body theory called the fluctuation 
exchange approximation (FLEX)lII is widely used in the paramagnetic phase. This is a self-consistent versioruof the 
ladder approximation, and has a nice feature of satisfying the conserving condition of Kadanoff and Baym.Ea One 
shortcoming of the FLEX is the lack of the vertex correction. For this reason, it may be interesting to extend the 
three-body scattering theory to the paramagnetic phase by setting up directly the integral equation instead of the 
Faddeev equation. 

We have also studied the instability of the Nagaoka state. The present theory is regarded as a natural extension of 
Roth's variational wave function, having improved the phase diagram of the Nagaoka state. This result may serve as 
a measure of the accuracy of the present theory. 

An accurate determination of the phase diagram is beyond the scope of this paper. If one uses more sophisticated 
variational wave functions, the area of the Najgaoka ground state can be reduced but remains finite. One of the best 
phase boundary was obtained by Wurth at al.fS who included more than 1000 terms with up to two particle- hole pair 
excitations in the variational calculation. Another interesting attempt was made by von der Linden and Edwards,E3 
who used the Slater determinant of the single-particle states of up-spin electrons infiuenced by the moving down-spin 
electron. The latter approach is regarded as an extension of the method of Richmond and Rickayzeuxj Both gave 
smaller areas of the Nagaoka state than the present result. On the other hand, Davis and FeldkampcJ pointed out 
that the method of Richmond and Rickayzen can be extended to calculate accurately the whole spectra of the single- 
particle excitations; one first diagonalize all ground and excited states in the presence of the static impurity of the 
down-spin electron and then diagonalize the hopping of the down-spin electron in this complete set of states. Of course 
the actual calculation has to be approximate by truncating the number of the excited states. Davis and FeldkampEj 
and one of the present authorsQ applied this method to finite-size systems of the one-dimensional Hubbard model, and 
calculated the whole spectral function. This type of calculations has not been attempted yet on the two-dimensional 
Hubbard model. 

So far, we have considered only the single-particle Green's function. The two-particle Green's function gives the 
poles of spin waves, whose |S£jftening may set another condition for the instability of the Nagaoka state. The phase 
diagram obtained by Okabecj along this line is not better than the present one. As a natural extension of the present 
theory to calculate spin waves, we need to consider four-body correlations. 

Finally we comment on the application to realistic situations of transition-metal compounds, in which five 3d- 
orbitals and the Hund-rule coupling are to be considered. The present theory is too complicated to treat such 
situations, and some simplifications have to be made. We have simplified the calculational scheme by using the local 
approximatioLL This type of calculations has been carried out for transition- metal mono-oxidesEil and La2Cu04 and 
Sr2Cu02Cl2,o which improves considerably the single-particle spectra of the band calculation based on the local 
density approximation and leads to good agreement with the photocmission experiments. 
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APPENDIX A: DERIVATION OF EQS. (2.40) AND (2.41) 



Expanding Di(k, z) and Z?2(k, z) in powers of 1/U, we rewrite Eqs. (2.21) and (2.22) as 



1 — n 1 

*Pp(Q, z; p) = -(1 -n) —{z + ep - /i) + -App(Q + p) 



^ I] |l + ^(^ + ep - /^) - -^[eq+p-q' + eq' - 2/^ - App(Q + p)] + • • •| 



*ph(Q, z- q) = 



X $ph(Q,z;q'), 

-[/7i-(z-eq + Ai) + iAph(Q-q) iA;,JQ-q,z) 



(Al) 



-1 



^(Q,^;q) 

1 A' 

I TTrF^ph 



J(Q,^;q)' 



1 A' (Q-q,z) 



J(Q,z;q) J(Q,2;q)2 



^S{- 



EQ+p'-q 



*Pp(Q,z;p'), 



where 



-^(Q, ^; q) = z - eq + /I Aph(Q - q), 



^pp(k) = T7 5I(^k-q' + eq' - 2//), 



q' 



^ph(k) = — ^(ek+p'- v)' 
p' 



(A2) 

(A3) 
(A4) 
(A5) 
(A6) 



Then we expand $pp(Q, z ; p) and $ph(Q, z; q) m powers of \/U as written down in Eqs. ( 2.3^ ) and ( p.3£ ). Substitute 
them into Eqs. ( |Ai| ) and (A2) and arranging the terms in each order of 1/C/, we obtain the following relations: 



Un 



U 



J(Q,z;q) J(Q,z;q)iV^ 
<I>W(Q,z;p) = -(l-n)-i-5]$W(Q,z;q' 



^E<^(Q'-;p'). 



$;,l)(Q,z;q) = -l 



N 
q' 

1 1 



(A7) 
(A8) 



J(Q,z;q)2n p^ 



A'{q-q,z) 



4h(Q-q.^)4E<HQ'^=p') 



n2j(Q,z;q)2 



p' 



^7(Q^^E'^^^(Q^^;p'), 



(A9) 



$W (Q, z; p) = ^{z + ep-fi) + -App(Q + p) 



1 1 

UN 



E 



£p - M (eq+p-q' + Cq' - 2^ - App(Q + p)) 



i>S(Q,^;q') 



lE<h'(Q,-;q')- 

q' 



(AlO) 



Note that p, p' are restricted inside the Fermi sphere and q, q' are outside the Fermi sphere. 

Equation ( |A8| ) indicates that ^pp (Q, z; p) i s in dependent of p. This fact and Eq.(A7) leads to <i>pp(Q,z;p) = —1. 
Substituting this relation into Eqs. (|A8[ ) and (A9), we have 



q' 



(All) 
(A12) 
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Here other terms in Eq. (A9) cancel out with each other. Equation ( A12| ) indicates that the q-dependence of 
*^ph (Q' "^i *l) comes only from J(Q, z; q). Therefore *&ph (Q, z] q) is written in a form 



*S(Q,^;q) 



c 



J(Q,z;q)' 
The constant C is determined by substituting Eq. (A13[) into Eq. (All|): 



N 2.^q J(Q,2;q) 



Substituting Eqs. (AlS) into Eq. ( [A12| ), we obtain 

p' 
Next, we rewrite Eq. ( |A10| ) as 

$«(Q, z- p) - -j^Y^^z - eQ+p_q. - eq. + ep + m)^S(Q, z- q') 



(A13) 



(AM) 



(A15) 



Un 



q' 



(A16) 



with the help of Eq. (All). Then we make sum with respect to p in both sides of Eq. (A16). The first term on the 
right hand side satisfies the relation 



]^ II(^ - eq+p-q' - eq' + ep + M)«'ph (Q: ^i ^) = riFiQ, z), 
P 

which is independent of q'. Using this relation, we have 



p' 



q' 



Substituting Eq. (A15) into this equation, we finally obtain 



q' 



U 



(A17) 



(A18) 



(A19) 



APPENDIX B: ADAPTION OF THE METHOD OF RICHMOND AND RICKAYZEN TO NUMERICAL 

CALCULATIONS 

Consider the situation that there exist N^ up-spin electrons and one down-spin electron. Neglecting the kinetic 
energy for the down-spin electron, we assume that it sits at origin. The up-spin electrons arc under the static repulsive 
potential localized at origin. Since there is no interaction working between up-spin electrons, we first calculate the 
single-particle energy by solving the potential problem, and then make N^ electrons occupy these energy levels 
according to the Pauli principle. 

Let the wave functions of single particles be 



lV'^"')=E"k Ik), n=l,2,...,7V, 



(Bl) 



where |k) represents the normalized state of plane wave with wave vector k. Then the energy eigenvalue i?*-"-* is given 
by 
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k' 

We solve numerically this equation by directly diagonalizing the corresponding matrix for systems of finite size (up 
to A^ = 60 X 60), although it is also possible to solve it by using the Green's function. Lowest Ne levels are occupied 
by up-spin electrons, and thereby the energy of the system is given by 

Etot = ^ +^1, (B3) 

n=l 

where Ei — 0, since we neglected the kinetic energy for down-spin electrons. The energy of the saturated ferromagnetic 
ground state for Ne + I up-spin electrons is given by E^ — J2lZi ' ^ki- If Etot is smaller than E^, the saturated 
ferromagnetic ground state is unstable. 



1 J. Hubbard, Proc. Royal Soc. London A276, 238 (1963). 
^ J. Kanamori, Prog. Theor. Phys. 30, 275 (1963). 
^ M. C. Gutzwiller, Phys. Rev. Lett. 10, 159 (1963). 

* J. A. Hertz and D. M. Edwards, J. Phys. F 3, 2174 (1973); D. M. Edwards and J. A. Hertz, J. Phys. F 3, 2191 (1973). 
^ H. Matsumoto, H. Umezawa, S. Seki, and M. Tachiki, Phys. Rev. B 17, 2276 (1978). 

^ J. Igarashi, J. Phys. Soc. Jpn. 52, 2827 (1983). 
■^ J. Igarashi, J. Phys. Soc. Jpn. 54, 260 (1985). 

* A. E. Ruckenstein and S. Schmitt-Rink, Int. J. Mod. Phys. B 3, 1809 (1989). 

The variational aspect of the three-body problem has discussed by E. G. Goryachev, Phys. Lett. A 166, 148 (1992). 
1" T. C. Hsu and B. Doucot, Phys. Rev. B 48, 2131 (1993). 
" Y. Nagaoka, Phys. Rev. 147, 392 (1966). 
^2 D. J. Thouless, Proc. Phys. Soc. 86, 893 (1965). 

1^ Y. Fang, A. E. Ruckenstein, E. Dagotto and S. Schmitt-Rink, Phys. Rev. B 40, 7406 (1989). 
" B. Doucot and X. G. Wen, Phys. Rev. B 40, 2719 (1989). 
^'^ A. Barbieri, J. A. Riera, and A. P. Young, Phys. Rev. B 41, 11697 (1990). 
^'^ W. von der Linden and D. M. Edwards, J. Phys.: Cond. Matter 3, 4917 (1991). 
" J. Hirsch, Phys. Rev. B 31, 4403 (1985). 

^** L. M. Roth, Phys. Chem. Solids 28, 1549 (1967); Phys. Rev. 186, 428 (1969). 
^^ S. R. Allan and D. M. Edwards, J. Phys. F 12, 1203 (1982). 

^" B. S. Shastry, H. R. Krishnamurthy and P. W. Anderson, Phys. Rev. B 41, 2375 (1990). 
^^ A. G. Basile and V. Elser, Phys. Rev. B 41, 4842 (1990). 

22 T. Okabe, Prog. Theor. Phys. 97, 21 (1997); Phys. Rev. B 57, 403 (1998). 

23 P. Richmond and R. Rickayzen, J. Phys. C 2, 528 (1969). 

2** T. Hanish, G. S. Uhlig and E. Miiller-Hartmann, Phys. Rev. B 56, 13960 (1997). 

25 L. D. Faddeev, Zh. Eksperim. Teor. Fiz. 39, 1459 (1960) [English transl.:Soviet Phys.-JETP 12, 1014 (1961)]. 

2^ This curve is very close to that of Roth's wave function. See also the curve "RES2" in Ref. 24. 

2^ N. E. Bickers and D. J. Scalapino, Annals of Physics, 193, 206 (1989). 

2** L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, Menlo Park, 1962). 

2^5 P. Wurth, G. S. Uhrig and Miiller-Hartmann, Ann. Phys. (Leipzig) 5, 148 (1996). 

3" L. C. Davis and L. A. Feldkamp, J. Appl. Phys. 50, 1944 (1979). 

3^ M. Takahashi and J. Igarashi, Ann. Phys. (Leipzig) 5, 247 (1996); Phys. Rev. B 54, 13566 ( 1996); Phys. Rev. B 56, 12818 

(1997). 
32 M. Takahashi and J. Igarashi, Phys. Rev. B 59, (1999) No.ll. 

FIG. 1. Feynman diagrams for the self-energy E|(Q, z) within the three-body scattering theory. Only three lines of Green's 
functions exist in the intermediate states. The broken lines represent the Coulomb interaction. 

FIG. 2. A sketch for a three-body state [Q; pq). 
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FIG. 3. Dispersion of spin-wave energy uiq along the line from q = (0,0) to {n,n) within the RPA. The density of up-spin 
electrons is n = 1, and the system size is A*' = 100 x 100. 

FIG. 4. Quasi-particle energy Eq and the lowest boundary of the continuum spectra as a function of Q. The calculation is 
carried out for a system of A = 50 x 50. The hole density and the Coulomb interaction are S{= 1 — n) — 0.2 and 4t/[/ = 0.05, 
respectively. 

FIG. 5. Quasi-particle energy Eq at Q = as a function of 4t/U. The hole density is 5 = 0.2 and the system size is 
N — 50 X 50; (a) the three-body scattering theory; (b) the ladder approximation. The line with lower energy in (b) is the 



contribution only from the first term of Eq. (2.2. 



FIG. 6. Quasi-particle energy Eq at Q = as a function of 5 in the U — oo limit. The system size is N = 100 x 100. The 
solid line represents the result of Roth's wave function (taken from Fig. 1 in Ref. 24). 

FIG. 7. Phase diagram for the square lattice. Open circles are for a system pf A = 40 x 40, and the solid squares are for a 
system of A = 50 x 50. The broken line is the curve "RES3" of Hanisch et alEj. The solid line with letters "SW" represents 
the boundary above which the instability of the RPA spin-waves become unstable. The solid line with letters "RR" represents 
the boundary given by the method of Richmond and Rickayzen. 
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